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Abstract. Low-frequency longitudinal oscillations of a flaring coronal loop are studied numerically. In the recent 
work of Nakariakov et al. A&A, 414, L25-L28 (2004) it has been shown that the time dependences of density and 
velocity in a flaring loop contain well-pronounced quasi-harmonic oscillations associated with a 2nd harmonic of 
a standing slow magnetoacoustic wave. In this work we investigate physical nature of these oscillations in greater 
detail, namely, we study their spectrum (using periodogram technique) and how does heat positioning affects the 
mode excitation. We found that excitation of such oscillations practically independent of the positioning of the 
heat deposition in the loop. Because of the change of the background temperature and density, the phase shift 
between the density and velocity perturbations is not exactly equal to the quarter of the period, it varies along 
the loop and is time dependent, especially in the case of one footpoint (asymmetric) heating. 
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1. Introduction 

Magnetohydrodynamic (MHD) coronal seismology is the 
main reason for studying waves in the solar corona. 
Also such studies are important in connection with 
coronal heating and solar wind acceleration problems. 
Observational evidence in the EUV coronal emission 
of coronal waves and oscillations is numerous (e.g. 
jOfman et al. 1999. Ofman &_Wang^J002J). Radio band 
observations also demonstrate various kinds of oscilla- 
tions (e.g., th e quasi-periodic pulsations, or QPP, see 
ijAschwanden lll987) for a review), usually with the pe- 
riods from a few seconds to tens of seconds. Also, decime- 
ter and microwave observations show much longer pe- 
riodiciti es, often in associa tion with a flare. For ex- 
ample, l)Wang &: Xie I l2000(l observed QPP with the 
periods of about 50 s at 1.42 and 2 GHz (in as- 
sociation with an M4.4 X-ray flare). Similar period- 
icities have been observed in the X-ray band (e.g., 
l)McKenzie fc Mullan I Il997l IXerekhov et al. Il2002|) ) and 
in the white-light emission associated with the stel- 
lar flaring loops ijMathioudakis et al. I l2003j) . A possi- 
ble interpretation of these medium period QPPs may 



be interpreted in terms of kink or torsional modes 
jZaitsev fc Stepanov 1989). In our previous, preliminary 
study jNakariako^ et al. .2004^) , we have outlined an alter- 
native, simpler, thus more aesthetically attractive, mecha- 
nism for generation of long-period QPPs. In this study we 
demonstrate in detail, that in a coronal loop an impul- 
sive (time-transient) energy release efficiently generates 
the second spatial harmonics of a slow magnetoacoustic 
mode. In particular, in the present study we study the 
spectrum (using periodogram technique) of these oscil- 
lations and how does heat positioning affects the mode 
excitation. 



Send offprint requests to: David Tsiklauri, 
e-mail: D . TsiklauriOsalf ord .ac.uk 



2. Numerical Results 

The model that we use to describe plasma dynami cs in a 
coronal loop is outlined in ijNakariakov et aril2004j) . Here 
we just add that, when solving numerically ID radiative 
hydrodynamic equations (infinite magnetic field approx- 
imation), and using a ID version of Lagrangian Re- map 
code (Arber et al. 2001) with the radiative loss limiters. 
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the radiative loss function was specified as 

T > 10^-6 
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which is ijRosner et al. I llQT^ law extended to a wider 
temperature range ijPeres et al. Ill982l IPriest lllQsJ) . 

We have used the same heating function as in 
ijNakariakov et al. Il2004j) . The choice of the temporal part 
of the heating function is such that at all times there is a 
small background heating present (either at footpoints or 
the loop apex) which ensures that in the absence of flare 
heating (when a, that determines flare heating amplitude, 
is zero) the average loop temperature stays at 1 MK. For 
easy comparison between that apex and footpoint heat- 
ing cases, with fix, Qp, flare heating amplitude at a given 
different value in each case (which ensures that with the 
flare heating on when a — 1 the average loop temperature 
peaks at about the observed value of 30 MK in both cases). 

In all our numerical runs, presented here, l/{2a'^) 
was fixed at 0.01 Mm-^^ which gives the heat deposition 
length scale, 0-^ = 7 Mm. Th is is a typical value deter - 
mined from the observations l)Aschwanden et al. l liool. 
The flare peak time was fixed in our numerical simula- 
tions at 2200 s. The duration of flare, at, was fixed at 
333 s. The time step of data visualization (which in fact 
is much larger than the actual time step, 0.034 s, in the 
numerical code) was chosen to be 0.5 s. 

2.1. Case of Apex (Symmetric) Heating 

In this subsection we complete analysis started in 
ijNakariakov et al. lEoO^ . namely for the same numeri- 
cal run we study spectrum of the detected oscillations at 
different spatial points. 

As it was pointed out in ( Nakariakov et al. 1 l2004ll . the 
most interesting fact is that we see clear quasi-periodic 
oscillations, especially in the second stage (peak of the 
flare) for time interval t = 2500 - 2800 s (cf. Fig. 1 in 
()Naka,ria,kov et a,l. 1 12004) '!. In fact, such oscillations are 
frequently seen during the solar flares o bserved in X-rays, 
8-20 keV, (e.g. iTerekhov et al. I tOO'j)) as well as stel- 
lar fla res observed in white-light fe.g. lMathioudakis et al. I 
(Hooi)). 

Before we start discussion of the physical nature of 
these oscillations, first let us recall for completeness what 
simple ID analytic theory of standing sound waves tells 
us: If we use ID, linearized, hydrodynamic equations with 
constant unperturbed (zero order) background variables. 
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Fig. 1. Case of apex (symmetric) heating: Periodogram 
(spectrum) of the velocity and density oscillatory compo- 
nent times series outputted in the following three points: 
loop apex (solid curve), 1/4 (dash-dotted curve) and 1/6 
(dashed curve) of the effective loop length (48 Mm), i.e. 
at s = 0,-12,-16 Mm. 
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Fig. 2. Case of apex (symmetric) heating: oscillatory com- 
ponents of time series outputted at ±6 Mm and ±18 Mm 
in the time interval of 2500-2800 s. The solid curve shows 
plasma number density in units of 10^^ cm-^. The dashed 
curve shows velocity normalized to 400 km s^^. 



the solutions for density, p and velocity, Vx, can be easily 
written as 



Tr / \ A I'mrCs \ . fniT 
i4;(s,t) = A COS ( — ^ — 1\ sm y'j^^ 



Apo . ( mTCs,\ friTT 



(1) 



(2) 



where Cg is a speed of sound, A is wave amplitude, L is 
loop length, n = 1,2,3, ... is a number of harmonic, and 
s is a distance along the loop. Note that (relative) phase 
shift between 14 and p is AP/P = -(7r/2)/(27r) = -1/4, 
where P is standing wave period. 
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In Fig. 1 we present periodogram (spectrum) of the 
velocity and density oscillatory component times series 
outputted in the following three points: loop apex, 1/4 
and 1/6 of the effective loop length (48 Mm), i.e. at 
s = 0,-12,-16 Mm. The choice first two points is ob- 
vious because we wanted to test whether simple analytic 
solution for ID standing sound waves (see below) is at all 
relevant in this case. The third point (1/6) was chosen 
arbitrarily (any spatial point along the loop where den- 
sity and velocity of the standing waves does not have a 
node would be equally acceptable). We gather from the 
graph that as expected for 2nd spatial harmonic of a 
standing sound wave in the velocity periodogram there 
are two clearly defined peaks and the largest peak corre- 
sponds to 1/4 of the effective loop length, while smaller 
peak corresponds to the 1/6. Note that at the loop apex 
periodogram gives (solid line is too close to zero to be 
seen in the plot). While density periodogram is exactly 
opposite to the velocity one. Largest peak corresponds to 
the loop apex, while 1/6 of the effective loop length cor- 
responds to a smaller peak, and 1/4 is quite small (close 
to zero). Location of the peaks are at about 0.0155 Hz i.e. 
period of the oscillation is 64 s. Period of a 2nd spatial 
harmonic of a standing sound wave should be 

P = L/Cs = L/{1.52x IO^Vt), (3) 

where T is plasma temperature measured in MK. If we 
su bstitute effecti ve loop length L — 48 Mm (see Fig. 2 
in ijNakariakov et al. 2004lj') ) and an a verage temperature 
of 25 MK (see top panel in Fig. 1 in ijNakariakov et ai~l 
l2004 h in the range of 2500-2800 s - quasi periodic os- 
cillations time interval we study) we obtain 63 s, which 
is quite close to the result of our numerical simulation. 
In fact, such a close coincidence is surprising bearing in 
mind that the theory does not take into account variation 
of background density and velocity in t ime, while we see 
from Fig. 1 in (^ Nakariakov et al'~ll2004l) that even within 
considered short interval during the flare 2500-2800 s all 
physical quantities vary in time significantly. To close our 
investigation of the physical nature of the oscillations, we 
study the phase shift between velocity and density oscil- 
lations, and compare our simulation results with the ID 
analytic theory. In Fig. 2 we plot oscillatory components of 
time series outputted at ±6 and ±18 Mm of plasma num- 
ber density in units of 10^^ cm~'^ and velocity, normalized 
to 400 km s^^ (note that sound speed at 25 MK is 760 km 
s^i which roughly is within the t = 2500-2800 s). These 
points where selected so that one symmetric (with respect 
to the apex) pair (±6 Mm) is close to the apex, while an- 
other pair (±18 Mm) is closer to footpoints. We gather 
from this graph: (A) clear quasi periodic oscillations are 
present; (B) that they are shifted with respect to each 
other in time; (C) in the case of the pair close to the apex 
(±6 Mm) the phase shift is quite close to the one predicted 
by the ID analytic theory; (D) in the case of the pair close 
to the footpoints (±18 Mm) the phase shift is somewhat 
different from the one predicted by the ID analytic theory. 
In the latter case the discrepancy can be attributed to the 



presence flows close to the footpoints. Main reason for the 
overall deviation is due to fact that analytic theory does 
not take into account variation of background density and 
velocity in time, and the fact that density gradients in 
the transition region are not providing perfect reflecting 
boundary conditions for the formation of standing sound 
waves. 

Yet another interesting result is that even with wide 
variation of the parameter space of the flare, its duration, 
peak average temperature, etc. we always obtained domi- 
nant 2nd spatial harmonic of a standing sound wave with 
some small admixture of 4th and sometimes 6th harmonic. 
Our initial guess was that this is due to the symmetric ex- 
citation of these oscillations (recall that we use apex heat 
deposition). In order to investigate this the issue of ex- 
citation further we decided to break symmetry and put 
heating source at one footpoint, hoping to see excitation 
of odd harmonics 1st, 3rd, etc. 

2.2. Case of Single Footpoint (Asymmetric) Heating 

In the case o f single footpoint heating we fix sq — 30 Mm 
in Eq.(l) in ijNakariakov et al'~ll2004j) . i.e. (spatial) peak 
of the heating is chosen to be at the bottom of the tran- 
sition region (top of chromosphere). Initially we run our 
numerical code without flare heating, i.e. we put a (in 
this manner we turn off flare heating). Eq = 0.02 erg cm^'^ 
s~^, was chosen such that in the steady (non-flaring) case 
the average loop temperature stays at about 1 MK. Then, 
we run our numerical code with flare heating, i.e. we put 
a = 1, and fix Qp at 1 x lO'^, so that it yields peak average 
temperature of about 30 MK. The results are presented 
in Fig. 3. We gather from the graph that during the flare 
apex temperature peaks at 38.38 MK while the number 
density at the apex at 3.11 x 10^^ cm~'^. In this case, as 
opposed to the case of symmetric (apex) heating, velocity 
dynamics is quite different. Namely, since the symmetry 
of heating is broken, there is non zero net flow through the 
apex at all times. However, as in the symmetric heating 
case, we again see quasi-periodic oscillations superimposed 
on the dynamics of all physical quantities (cf time interval 
of < = 2400—2700 in Fig. 3). In what follows we do analysis 
of their physical nature as in the previous subsection. 

In Fig. 4 we present time-distance plots of velocity and 
density in the time interval 2400-2700s, where the quasi- 
periodic oscillations are most clearly seen. Here we again 
subtracted the slow varying with respect to oscillation pe- 
riod background. The picture is quite different from the 
cas e of apex (symmetric) h eating (compare it with Fig. 2 
in llNakariakov et al. I2OO4I) ). It is looks more complex be- 
cause in this case the node in the velocity (at the apex) 
moves back and forth periodically along the apex, and 
for larger times {t > 2550 s) stronger flows (darker bands 
with at a slope) are now present. However, physical nature 
of the oscillations remains mainly the same - 2nd spatial 
harmonic of a sound wave (though with oscillating node 
along the apex). 
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To investigate this further we plot in Fig. 5 peri- 
odogram (spectrum) of the velocity and density oscilla- 
tory component times series outputted in the following 
three points: loop apex, 1/6 and 1/4 of the effective loop 
length. We gather from the graph that the periodogram 
(spectrum) is more complex than in the case of apex (sym- 
metric) heating. Namely, in the velocity periodogram at 
the apex we see a peak with frequency of the oscillation 
different (higher) than that of 2nd spatial harmonic of a 
standing sound wave. Actually that is the frequency with 
which node of the velocity oscillates (see discussion in the 
previous paragraph) . It has nothing to do with the stand- 
ing mode, and it is rather dictated by the excitation condi- 
tions in the dynamics resonator, which the considered loop 
obviously is. Let us analyze now how this periodogram 
compares with the ID analytic theory. Peak in the peri- 
odogram corresponding to 1 /6 of the effective loop length 
(dashed hne) corresponds to frequency of about 0.017 Hz, 
i.e. period of the oscillation is 59 s. Again, the period of 
a 2nd spatial harmonic of a standing sound wave should 
be calculated from Eq. (3). If we substitute effective loop 
length L = 48 Mm (see Fig. 4) and an average temper- 
ature of 26 MK (see top panel in Fig. 3 in the range of 
2400-2700 s - quasi periodic oscillations time interval we 
study) we obtain 62 s, which is quite close to the result of 
our numerical simulation (59 s). 

Next, we study the phase shift between velocity and 
density oscillations, and compare our simulation results 
with the ID analytic theory. In Fig. 6 we produce a plot 
similar to Fig. 2, but for the case of asymmetric heat- 
ing. The deviation which is greater than in the case apex 
(symmetric) heating can again be attributed to the over- 
simplification of the ID analytic theory, which does not 
take into account time variation of the background phys- 
ical quantities and imperfection of the reflecting bound- 
ary conditions (see above). And, more importantly, in the 
asymmetric case strong flows are present throughout of 
the flare simulation time, thus, if, say, linear time depen- 
dence would be assumed, which is relatively good within 
considered short interval during the flare 2400-2700 s, the 
Eqs(l)-(2) would be modified such that phase shift would 
vary secularly in time, which is similar to what we see in 
Fig. 6. 

Yet another interesting observation comes from the fol- 
lowing argument: in a steady ID case analytic theory pre- 
dicts that the phase shift between the density and velocity 
should be (A) zero for for a propagating acoustic wave and 
(B) quarter of a period for the standing acoustic wave. 
Since in the asymmetric case strong flows are present, we 
see less phase shift between the velocity and density in 
Fig. 6 as one would expect. 

Thus, the results of th e present study provide further 
more definitive than in ( Nakariakov et al~l l2004l) proof 
that these oscillations indeed are 2nd spatial harmonic of a 
standing sound wave. However, present work also revealed 
that in the case of single footpoint (asymmetric) heating 
the physical nature of the oscillations is more complex as 
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Fig. 3. Case of single footpoint (asymmetric) heating: 
Average temperature, temperature at apex, and number 
density at apex as a function of time. 
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Fig. 4. Time-distance plots of the velocity and density 
oscillatory components in the time interval of 2400-2700s 
for the case of single footpoint (asymmetric) heating. 



the node in the velocity oscillates along the apex and net 
flows are present. 



3. Conclusions 

Initially we have used ID radiative hydrodynamics loop 
model which incorporates the effects of gravitational strat- 
ification, heat conduction, radiative losses, added ex- 
ternal heat input, presence of Helium, hydrodynamic 
non-lin earity, and bulk Brag inskii viscosity to simulate 
flares (jTsiklauri et al. 1 120041) . As a byproduct of that 
study, in practically all our numerical runs we have de- 
tected quasi periodic osc illations in all physical quantities 
((Nakariakov et al'~ll2004l) . In fact, such oscillations are fre- 
quently seen during the solar flares observed in X-rays, 
8-20 keV (e.g. i^rckhov et al. (2 003)) as weU as stel- 
lar fla res observed in white-light (e.g. lMathioudakis et al. I 
((2003!) ). Our present analysis shows, in detail, that quasi 
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Fig. 5. As is Fig. 1 but for the case of single footpoint 
(asymmetric) heating. Time interval here is 2400-2700s. 
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Fig. 6. As in Fig. 2, but for the case of single footpoint 
(asymmetric) heating. Time interval here is 2400-2700s. 

periodic oscillations seen in our numerical simulations 
bear r nany simila r features as the observed ones. In sum- 
mary ijNakariakov et al. 11200^ and the present study es- 
tablished the following features: 

— We show that the time dependences of density and 
temperature in a flaring loop contain well-pronounced 
quasi-harmonic oscillations associated with standing 
slow magnetoacoustic modes of the loop. 

— For a wide range of physical parameters, the dominant 
mode is the second spatial harmonic, with a velocity 
oscillation node and the density oscillation maximum 
at the loop apex. This result is practically independent 
of the positioning of the heat deposition in the loop. 

— Because of the change of the background temperature 
and density, and the fact that density gradients in 
the transition region are not providing perfect reflect- 
ing boundary conditions for the formation of standing 
sound waves, the phase shift between the density and 
velocity perturbations is not exactly equal to the quar- 
ter of the period. 



— We conclude that the oscillations in the white light, the 
radio and X-ray light curves observed during solar and 
stellar flares may be produced by the slow standing 
mode, with the period determined by the loop tem- 
perature and length. 

— For a typical solar flaring loop the period of oscillations 
is shown to be about a few minutes, while amplitudes 
are typically few percent. 

The novelties brought about by this study are that 
by studying the spectrum and phase shift of these oscil- 
lations we provide more definite proof that these oscilla- 
tions are indeed 2nd harmonic of a standing sound wave, 
and that the single footpoint (asymmetric) heat position- 
ing still produces 2nd spatial harmonic, though it is more 
complex than the apex (symmetric) heating (due to the 
presence of flows). 
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